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In this Analysis we investigated the possibility of the measurement of Higgs self-coupling at ILC 
through the process e + + —¥ ZHH using fast simulation data. So far two combinations of decay 
modes: Z —¥ qq, H — > bb, H — > WW* and Z —¥ ll,H —¥ bb,H —¥ bb were studied. Our preliminary 
results show that it is very challenging to suppress the huge standard model backgrounds effectively. 



I. INTRODUCTION 



It is well accepted that the discovery of a Higgs-like boson is not enough to fully understand the mechanism of 
electro-weak symmetry breaking (EWSB) and mass generation. The Higgs self-coupling can be a non-trivial probe of 
the Higgs potential and probably the most decisive test of the EWSB mechanism. In the standard model framework, 
the Higgs potential V(&) = \(& 2 — ^v 2 ) 2 , where $ is an isodoublet scalar field and v rj 246 GeV is the vacuum 
expectation value of its neutral component, is uniquely determined by the self-coupling A. Obviously, determination 
of the Higgs mass, which satisfies m 2 H = 2\v 2 at tree level, can provide an indirect information about the self- 
coupling. The measurement of the trilinear self-coupling Xhhh = 6Aw offers direct independent determination of the 
Higgs potential shape, which is the topic of this analysis. 

The trilinear Higgs self-coupling can be measured at ILC through the two leading processes: double Higgs-strahlung 
and WW fusion 043: which are shown in FigJT] The former is expected to dominate around the center of 
mass energy of 500 GeV and the latter to take it over at higher energy. In this analysis we focus on the double 



Higgs-strahlung process e" 1 



ZHH for the Higgs mass of Mh = 120 GeV and the center of mass energy of 



^fs = 500 GeV with the integrated luminosity 2 ab 

Depending on the different decay modes of Z and H, there are different methods to identify the signal events. Table 
I shows several most promising combinations of decay modes for e + + e~ — > ZHH and their branching ratios. Modes 
1 and 3 are studied in Ref. [8|. We study the other two modes in this analysis. 
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FIG. 1: Leading processes involving trilinear Higgs self-coupling: (Left) Double Higgs-strahlung; (Right) WW fusion. 



TABLE I: Most promising modes for e+ + e~ ->• ZHH 



Decay Mode 


Z -> 


Hi ->■ 


H 2 ^ 


Branching Ratio 


1 


qq 


bb 


bb 


34% 


2 


qq 


bb 


WW* 


14% 


3 


vv 


bb 


bb 


9.8% 


4 


11 


bb 


bb 


4.9% 
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II. SIMULATION 



The simulations of signal events (e+ + e~ — > ZHH) and possible background events (e + + e~ — > 
tt, ZZZ, W + W~Z, ZZ, ZH) were done by Physsim In Physsim the helicity amplitudes arc calculated by the 
HELAS library [10| , The phase space integration and the four momenta generation are performed by BASES/SPRING 
plj . Parton showering and hadronization are carried out by PYTHIA6.4 [12], where final-state t leptons are decayed 
by TAUOLA 13| in order to handle their polarizations correctly. The detector simulation was done by JSFQuick- 
Simulator, which implements the GLD geometry and other detector-performance related parameters [14| . 

It is worth mention of that the simulations were performed without the beam polarization but with the initial-state 
radiation, beam width and bcamstrahlung. Then the cross sections used here are shown in Tabic II. An integrated 
luminosity of 2 ab -1 is assumed in this analysis. 



TABLE II: Cross sections of the related processes 



Process 


e+ + e~ -> ZHH 


e+ + e~ ->• tt 


e+ + e~ -)• ZZZ 


e+ +e~ -> W+W-Z 


e+ + e- ^ ZZ 


e+ + e" -> Zif 


Cross section 


152 ab 


530 fb 


800 ab 


36 fb 


515 fb 


70 fb 



III. ANALYSIS 



A. e + +e~ ^ ZHH (qq){bb)(WW*) 



The full hadronic decays of W and VF* were investigated. In this mode the final state of a candidate signal event 
contains of 8 jets, two of which are b jets. To select the signal events, first we find all the good tracks and require the 
number of tracks be greater than 20. We then try to combine tracks with a small Y value to a current jet cluster, 
where the Y value between two momenta pi,P2 is defined as Y{p\,p2) = M J? 1 ' ; with M(pi,p2) being the invariant 
mass of pi,P2 and E vis the total visible energy. We continue the jet clustering until there are 7 jets left, because the 
two jets coming from W* are very close to each other which means the Y- value between them is very small, thereby 
being likely to be clustered as one jet. At this point we calculate the Y values for all the pairs from these 7 jets and 
choose the minimum denoted by Y cut . The Y cut distributions of signal events and background events (here we consider 
the tt events as background) are shown in Figj2] The 7 jets are combined by minimizing the \ 2 which is defined as 
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(M(b,b)-M H ) 2 (M(W, W*) — Mjj) 2 (M(q,q)-M z ) 2 (M(q,q') - M w f 



where M{q,q') is the reconstructed invariant mass of jet q and jet q 1 , Mjj, Mz and Myy are the mass of H, Z and 
W, respectively, and <thi , &h 2 7 ozandaw are their corresponding mass resolutions. 

In order to further suppress the background, we require that x 2 < 20, 90GeV < M(H\) < 130GeV, HOGcV < 
M(H 2 ) < 150GeV, 70GcV < M(Z) < 110GeV,y cut > 0.0076, where the asymmetry of two Higgs mass requirement 
is due to their different decay modes. The preliminary result of this cut-based analysis is shown in Table III. Though 
we can still add other cuts like b tagging requirement, the signal events will become too few to be observed. It seems 
very challenging to reject the huge tt background in this mode. 

We are going to investigate the semi-lepton decays of W and W* . 



TABLE III: Cut statistics of e+ + e~ -> ZHH -)■ (qq)(bb)(WW*) 
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Process 


ZHH -> (qq)(bb)(WW*) 


tt 


theoretical 


18.3 


1062000 


pre-selection 


12.6 


483949 


Y 1 < 20 


5.2 


65144 


90GeV < M Hl < 130GW 


5.1 


63157 


HOGeV < M H , < ISOGeV 


3.6 


36670 


90GeV < M z < UOGeV 


3.5 


34359 


Y cut > 0.005 


2.3 


8454 


Y cut > 0.0076 


1.1 


2644 




FIG. 2: Distribution of Y C ut, where black is for signal and red is for tt background. 



B. e+ + e~ -> ZHH ->■ (ll)(bb)(bb) 



In this mode a candidate signal event contains two leptons and four b jets, where we only consider the Z boson 
decaying into e + e~ and The two isolated charged lepton tracks are required to have an energy greater than 

20 GeV and the energy deposited in the cone of 20° around each lepton track be less than 20 GeV. We then force the 
other tracks to four jets and combine the four jets by minimizing the \ 2 defined by 

2 _ (M(b, b) - M H ) 2 (M(b, b) - M H ) 2 (M(l, Q - M z ) 2 

Table IV shows that 15.4 signal events survived the pre-selection but with thousands times more background events 
left. In order to reject the background effectively, while keeping a reasonable signal efficiency, we used the neural net 
method MLP in the TMVA package which gives some useful classifiers. Here we mainly consider the five kinds 
of backgrounds that are shown in Table IV. First we separately do the neural net analysis between the signal and 
each of the five kinds of backgrounds. For each event we can get five classifiers to separate signal and backgrounds. 
Figure |3] histograms the classifiers obtained by the MLP method for the signal and it samples. 



We then add some more cuts on the five classifiers denoted by mva_tt, mva_zzz, mvajwwz, mva_zz and mvajzh as 
shown in Table IV. We further impose cuts on the Z mass, Y cut , and require b tagging, which is based on the number 
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TMVA response for classifier: MLP 
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FIG. 3: The classifier obtained by neural net training for signal and tt background. 



of tracks with 2.5a separation from the interaction point. Our preliminary result is listed in Table IV. The final cut 
is applied with the neural net for the signal and the ZZZ background after all of the above cuts. We end up with 
3 signal events with 0.82 ZZZ events left, while four backgrounds are eliminated. The result shows that the ZHH 
events can be observed in this mode with the significance jSrm ~ 1.5a. 



TABLE IV: Cut statistics of e+ + e~ -> ZHH -)■ (ll)(bb){bb) 



Process 


ZHH 


tt 


ZZZ 


WWZ 


ZZ 


ZH 


generated 


1M 


4.5M 


500K 


750K 


1.25M 


250K 


theoretical 


304 


1062000 


1600 


72300 


1030000 


140000 


pre-selection 


15.4 


9023 


125 


1943 


3560 


1618 


mvaJt > 0.98 














mvajwwz > 1.0 














mva,-zz > 0.97 


11.7 


312 


12.9 


12.7 


16.5 


5.6 


mva_zh > 0.97 














mva-zzz > 














70GeV < Mz < HOGeV 


9.7 


106 


11.7 


7.5 


16.5 


0.56 


Y cut > 0.015 


9.1 


91.3 


10.6 


6.9 


6.6 





2b(H 1 ){N off >0) 


6.3 


28 


5.5 


1.8 








2b(H 2 )(N off > 1) 


3.5 


0.71 


2.3 











mvajzzz > 0.86 


3.0 





0.82 












IV. SUMMARY 



The two modes, e+ + e" -> ZHH (qq){bb){WW*) and e+ + e" ZHH (ll)(bb)(bb), were investigated for 
the purpose of the measurement of the trilinear Higgs self-coupling at ILC for Mh = 120 GeV, y/s = 500 GeV and 
the integrated luminosity of 2 ab" 1 . The former mode is very difficult to use for the signal observation, while the 
latter mode can be useful to observe the self-coupling. 



5 



Acknowledgments 

We would like to thanks all the members of the ILC physics subgroup [Tg| for useful discussions. This study 
is supported in part by KEK, Center of High Energy Physics, Tsinghua University and the JSPS Core University 
Program. 



[1] G. Gounaris, D. Schildknecht and F. Renard, Phys. Lett. B83 (1979) 191 and (E) 89B (1980) 437; V. Barger, T. Han 
and R.J.N. Phillips, Phys. Rev. D38 (1988) 2766. 

V.A. Ilyin, A.E. Pukhov, Y. Kurihara, Y. Shimizu and T. Kaneko, Phys. Rev. D54 (1996) 6717. 
A. Djouadi, W. Kilian, M. Muhlleitner and P.M. Zerwas, Eur. Phys. J. C10:27-43, 1999. 
F. Boudjema and E. Chopin, Z. Phys. C73 (1996) 85. 
V. Barger and T. Han, Mod. Phys. Lett. A5 (1990) 667. 
A. Dobrovolskaya and V. Novikov, Z. Phys. C52 (1991) 427. 

D.A. Dicus, K.J. Kallianpur and S.S.D. Willenbrock, Phys, Lett. B200 (1988) 187; A. Abbasabadi, W.W. Repko, D.A. 
Dicus and R. Vega, Phys. Re v. D38 (1988) 2770; Phys. Lett. B213 (1988) 386. 
Y. Takubo, arXiv:0907.0524vl. 



http: / / acfahep.kek.jp / subg/ sim/softs.html 



H. Murayama, I. Watanabe, K. Hagiwara, KEK-91-11, (1992) 184. 
T. Ishikawa, T. Kaneko, K. Kato, S. Kawabata, Comp. Phys. Comm. 41 (1986) 127. 
T. Sj ostrand, Comp. Phys. Comm. 82 (1994) 7 4. 
|http: / /wasm.home.cernch /wasm / goodies.html 



GLD Detector Outline Document, arXiv:physics/0607154 



A. Hoecker, P. Speckmayer, J. Stelzer, J. Therhaag, E. von Toerne, H. Voss et al, arXiv:physics/0703039 ^5. 
|http: / /www-jlc.ke k.jp/subg/physics/ilcphys 



